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Homogenization for advection-diffusion in a 
perforated domain 

p. H. Hayncs^, V. H. Hoang^, J. R. Norris^ and K. C. Zygalakisf 



Abstract 

The volume of a Wiener sausage constructed from a diffusion process 
with periodic, mean-zero, divergence-free velocity field, in dimension 3 
or more, is shown to have a non-random and positive asymptotic rate of 
growth. This is used to establish the existence of a homogenized limit for 
such a diffusion when subject to Dirichlet conditions on the boundaries 
of a sparse and independent array of obstacles. There is a constant ef- 
fective long-time loss rate at the obstacles. The dependence of this rate 
on the form and intensity of the obstacles and on the velocity field is in- 
vestigated. A Monte Carlo algorithm for the computation of the volume 
growth rate of the sausage is introduced and some numerical results are 
presented for the Taylor-Green velocity field. 
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1 Introduction 

Wc consider the problem of the existence and characterization of a homo- 
genized limit for advection-difFusion in a perforated domain. This prob- 
lem was initially motivated for us as a model for the transport of water 
vapour in the atmosphere, subject to molecular diffusion and turbulent 
advection, where the vapour is also lost by condensation on suspended 
ice crystals. It is of interest to determine the long-time rate of loss and 
in particular whether this is strongly affected by the advection. In this 
article we address a simple version of this set-up, where the advection is 
periodic in space and constant in time and where the ice crystals remain 
fixed in space. 

Let K he a, compact subset of M.'^ of positive Newtonian capacity. We 
assume throughout that d ^ 3. Let p G (0, oo). We consider eventually 
the limit p 0. Construct a random perforated domain D C M.'^ by 
removing all the sets K + p, where p runs over the support P of a 
Poisson random measure p, on M'' of intensity p. Let u be a Z''-periodic, 
Lipschitz, mean-zero, divergence-free vector field on W^. Our aim is to 
determine the long-time behaviour, over times of order = p~^, of 
advection-diffusion in the domain D corresponding to the operatoi0 

C = ^A + v{x).\7 

with Dirichlet boundary conditions. It is well known (see Section 
that the long-time behaviour of advection-diffusion in the whole space 
R'' can be approximated by classical, homogeneous, heat-flow, with a 
constant diffusivity matrix a = a{v). The effect of placing Dirichlet 
boundary conditions on the sets K + p is to induce a loss of heat. The 
homogenization problem in a perforated domain has been considered 
already in the case of Brownian motion , : [13 : 15 1 and Brownian 



motion with constant drift Q . The novelty here is to explore the possible 
interaction between inhomogeneity in the drift and in the domain. We 
will show that as p — )• there exists an effective constant loss rate 
X{v, K) in the time-scale . We will also identify the limiting values of 
r'^~'^X{v,rK) as r — > and ?' — )■ oo and we will compute numerically 
this function of r for one choice of v and K. 

Fix a function / £ L'^{M.'^). Write u = u{t,x) for the solution to 
the Cauchy problem for C in [0,cx)) x D with initial data /, and with 

^ All results to follow extend to the case of the operator ^ divaV -I- ii(a;).V, where 
a is a constant positive-definite symmetric matrix, by a straightforward scaling 
transformation. We simplify the presentation by taking a = I. Results for the 
case a = e^I are stated in Section [7] for easy reference. 
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Dirichlet conditions on the boundary of D. Thus, for suitably regular K 
and /, u is continuous on [0, oo) x D and on (0, oo) x D, and is C^'^ on 
(0, oo) X D] we have u(0, x) = f{x) for aU a; G D and 

— ~ ^Aw + v{x).Vu on (0, oo) x D. 

We shall study the behaviour of u over large scales in the limit p — ?> 0. 
Our analysis will rest on the following probabilistic representation of u. 
Let X be a diffusion process in R'', independent of /x with generator C 
starting from x. Such a process can be realised by solving the stochastic 
differential equation 

dXt = dWt + v{Xt) dt, Xq = x (1.1) 

driven by a Brownian motion W in R''. Set 

T = M{t^O:Xt eK + P}. 

Then 

u{t,x)=E, (/(XOl{T>t}|A*)- 

The key step is to express the right hand side of this identity in terms 
of an analogue for X of the Wiener sausage. Associate to each path 
7 € C([0,oo),M'^) and to each interval / C [0,oo) a set 5'f (7) C R'^ 
formed of the translates of K by 74 as t ranges over /. Thus 



5f (7) = Utei{K + 7t) = {x e R'^ : x - 7f e K for some t e I}. 

'(at]( 



Write Sf^ for the random set Sf^^JX) and write \S^\ for the Lebesgue 



volume of 5/^. We call Sj^ the diffusion sausage or {X, K)-sausage and 
refer to K as the cross section. Then T > t if and only if fi{S^) = 0, 
where K ~ {~x : x G K}. Hence 

u{t,x) ==E^ (/(Xt)l^^(gA-)^oj m) 
and so, by Fubini, we obtain the formulae 

E{u{t, x)) = (/(Xi) cxp(-p|5f I)) (1.2) 

and 

E{u{t, xf ) = E., [f{Xt)f{Yt) exp(-p|5f (X) U 5f (r)|)) (1.3) 

where Y is an independent copy of X. 

In the next section we review the homogenization theory for C in the 
whole space. Then, in Section[3]we show, as a straightforward application 
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of Kingman's subadditive ergodic theorem, that the sausage volume \S^\ 
has almost surely an asymptotic growth rate 7(1;, iiT), which is non- 
random. In Section 3] we make some further preparatory estimates on 
diffusion sausages. Then in Section [S] we identify the limiting values of 
r'^~'^^{v, rK) as r — and as r — >■ 00. In SectionlHl we use the formulae 
(|1.2p . (|1.3p to deduce the existence of a homogenized scaling limit for 
the function u, and we prove a corresponding weak limit for the diffusion 
process X and the hitting time T. We shall see in particular that for large 
obstacles it is the effective diffusivity a which accounts for the loss of heat 
in the obstacles. On the other hand, when the obstacles are small, the 
loss of heat is controlled instead by the molecular diffusivity, even over 
scales where the diffusive motion itself is close to its homogenized limit. 
Some results for non-unit molecular diffusivity are recorded in Section 
[71 Finally, in Section |51 we describe a new Monte Carlo algorithm to 
compute the volume growth rate for the (X, iir)-sausage, and hence the 
effective long-time rate of loss of heat. We present some numerical results 
obtained using the algorithm which interpolate between our theoretical 
predictions for large and small obstacles. 



2 Review of homogenization for diffusion with 
periodic drift 

There is a well known homogenization theory for /^-diffusion in the whole 
space R'^. See [H, 0j HI- We review here a few basic facts which 



provide the background for our treatment of the case of a perforated do- 
main. Our hypotheses on v ensure the existence of a periodic, Lipschitz, 
antisymmetric 2-tcnsor field (5 on such that ^ div/3 = v. So we can 
write C in the form 

£ = idiv(/ + /3(.T))V. 

Then C has a continuous heat kernel p : (0, and there 

exists a constant C < 00, depending only on the Lipschitz constant of 
V, such that, for all i, x and y, 

C-^cM~C\x-v\''/t}^p{t,x,y) ^ Cexp{-|x-y|VC^t}- (2.1) 



Moreover, C may be chosen so that there also holds the following Gaus- 
sian tail estimate for the diffusion process X with generator C starting 
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from x: for alH > and (5 > 0, 

sup - > (5 ) 5$ Ce-^'/^*. (2.2) 



The preceding two estimates show a quahtative equivalence between 
X and Brownian motion, vahd on all scales. On large scales this can be 
refined in quantitative terms. Consider the quadratic form q on given 

by 



g(0 = inf / 1$ - divx + /3V6l|^ dx 

where the infimum is taken over all Lipschitz functions 9 and all Lipschitz 
antisymmetric 2-tensor fields x on the torus T'' = /U^. The infimum 
is achieved, so there is a positive-definite symmetric matrix a such that 

g(0 = (e,a"'0- 

The choice 6* = and x = shows that a ^ /. As the velocity field 
V is scaled up, typically it is found that a also becomes large. See for 
example for further discussion of this phenomenon. 

We state first a deterministic homogenization result. Let / e L'^{W^) 
and (J S (0, oo) be given. Denote by u the solution to the Cauchy problem 
for C in R'' with initial data /(./cr) and set u^'^\t, x) ~ u{a^t, ax). Then 

\u'^''\t,x)-u{t,x)\^dx^Q (2.3) 

as cr — > oo, for all i ^ 0, where u is the solution to the Cauchy problem 
for i div aV in R'' with initial data /. 

In probabilistic terms, we may fix x G R'' and a E (0, oo) and consider 
the £-diffusion process X starting from ax. Set X^"^ = (T-^X^2t. Then 



it is known [13[ that 

X'-"'^ X, weakly on C([0, oo), R'') (2.4) 

where X is a Brownian motion in R"^ with difFusivity a starting from 
X. The two homogenization statements are essentially equivalent given 
the regularity implicit in the above qualitative estimates, the Markov 
property, and the identity 

u^'^\t,x)=E{f{X^/^)). 
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3 Existence of a volume growth rate for a diffusion 
sausage with periodic drift 

Recall that the drift v is Z'^-periodic and divergence-free. 

Theorem 3.1 There exists a constant 7 — "f{v, K) G (0, 00) such that, 
for all X, 

lim — - — = 7, fx-almost surely. 
t— >oo t 

Proof Write tt for the projection T''. Since v is periodic, the 

projected process n{X) is a diffusion on T''. As v is divergence-free, the 
unique invariant distribution for 7r(X) on T'' is the uniform distribution. 
The lower bound in (|2.ip shows that the transition density of tt{Xi) on 
T is uniformly positive. By a standard argument 7r(X) is therefore uni- 
formly and geometrically ergodic. Consider the case where Xq is chosen 
randomly, and independently of W , such that 7r(Xo) is uniformly dis- 
tributed on T''. Then 7r(X) is stationary. For integers ^ m < n, define 
Vm,n = l'5'(^^,„]l- Then Vi^n ^ Vi,^ -I- Kn,„ whenever ^ / < m < n. 
Since Lebcsgue measure is translation invariant and 7r(X) is stationary, 
the distribution of the array {Vm+k.n+k : ^ m < n) is the same for 
all fc 5^ 0. Moreover Vm,n is integrable for all m, n by standard diffusion 
estimates. Hence by the subadditive ergodic theorem Q we can conclude 
that, for some constant 7 5^ 0, 

lim J — — = 7, almost surely. 

n— >oo Ti 

The positivity of 7 follows from the positivity of cap(_R') using Theorem 
15. II below. 

Let Pj; be the probability measure on C([0,oo),R'') which is the law 
of the process X starting from x. Set 

g{x) ^fj lim ^ = 7V .9(2^) = 1™ = 7 I • 



Then g is periodic and g ~ g. We have shown that 

g{x) dx ~ 1. 



Hence g{x) = 1 for Lebesgue-almost-all x. But then by the Markov 
property, for every x, 

g{x)=g{x)^l p{l,x,y)g{y)dy = 1 
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which is the desired almost-sure convergence for discrete parameter n. 
An obvious monotonicity argument extends this to the continuous para- 
meter t. □ 



4 Estimates for the diffusion sausage 

We prepare some estimates on the diffusion sausage which wih be needed 
later. These are of a type well known for Brownian motion Q and extend 
in a straightforward way using the qualitative Gaussian bounds (|2.ip and 

Lemma 4.1 For all p e [l,oo) there is a constant C(j),v,K) < oo 
such that, for all t ^ and all x € R'', 

Proof Reduce to the case t = 1 by subadditivity of volume and L^- 
norms and by the Markov property. The estimate then follows from 
((2?2]) since 

|5f I s; LJd f rad(i^) + sup \Xt - x\ 

□ 

Lemma 4.2 There is a constant C{v,K) < oo with the following 
property. Let X and Y be independent C-diffusions starting from x. For 
all t ^ 1 and all x G K'', for all a,b ^ 0, 

(^('^aa+i)t](^) n ^ 0) ^ C{a + (4.1) 

and, when b ^ a + 1, 

P. {S(i,^ia+i)t] (X) n (X) ^ 0) < C(6 - a - l)-'^/^. (4.2) 

Proof Wc write the proof for the case t ~ 1. The same argument applies 
generally. There is alternatively a reduction to the case i = 1 by scaling. 
Assume that b ^ a + 1. Write J-t for the cr-algcbra generated by {Xg : 
^ s ^ t) and set 

i?Q = sup \Xt-Xa+i\, Rb^ sup \Xt-Xb\, Z = Xb-Xa+l- 
Then, by (1^ . 

P,(i?6 ^ \Z\/i\Tb) ^ Ce-l^l'/9C^ 
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so, using (|2.ip . 

P,(i?b ^ \Z\/3) CE,(e-l^l'/9^) C(5 - a - 
On the other hand, by agam, 

P,(i?a ^ |^|/3| J-a+i) C(6 - a - if^^Rt 
so, using (j2.2p . 

P,(i?a ^ |Z|/3) «;C(6-a-l)'^/^ 

Moreover (|2.ip gives also 

P:r(2rad(/v) ^ |Z|/3) =^ C{b~a^lf^^. 

Now if n ^ then either i?, > |Z|/3 or i?^ ^ 

|Z|/3 or 2rad(iir) ^ \Z\/3. Hence the preceding estimates imply (|4.2I) . 
The proof of (|4.ip is similar, resting on the fact that — Y}, has density 
bounded by C(a + 6)^''/^, and is left to the reader. 

Lemma 4.3 yls t — > oo, we have 



□ 



sup 







and 



|5f(x)n5f(r)| 



t 



sup Ea- 

a: 

where Y is an independent copy of X 
Proof Note that 



^ 



|5'^^t+i]| < \Sl^i \ ^ |5'f I + |S'(^ t^i]| 



Given Lemma l4.ll the first assertion will follow if we can show that, as 
t — oo, 



sup Ej. 



■7 



0. 



But by the Markov property and using (|2.ip , 



E, 



^ C 



E„ 



[0,1]'' 



\S, 



K\ 



p(l,x,y)Ey 
dy^Q 



dy 
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as i — > oo, where we used the almost-sure convergence \S^\/t — >■ 7 when 
7r(Xo) is uniform, together with uniform integrability from Lemma l4.ll 
to get the final limit. 

For the second assertion, choose q G (1,3/2) and p E (3,03) with 
1/p+l/q = 1. Then, for j,k ^ 0, by Lemmas 14.11 and 14.21 there is a 
constant C(p, v, K) < 00 such that 

E,(|5(^^^.+i](x)n5(^,,+i](y)|) 

C{j + k)"^/'^". 
So, as n — >■ cxj, we have 

•|5,f(x)n5f(r)| 



1/9 



E, 



n 



, E.(|^f(x)|) ^ -^-^ E.(|gg,^,j(x)n5(l,^,j(y)|) 

j=i k=o 



□ 



5 Asymptotics of the growth rate for small and 
large cross-sections 

We investigate the behaviour of the asymptotic growth rate "f{v, rK) of 
the volume of the {X, ri4')-sausage S^^ in the limits r — > and r — >■ 00. 
Recall the stochastic differential equation (jl.ip for X and recall the 
rescaled process X^"") from Section[21 Set W^"^ = a^^^W^^f Then W'^"^ 
is also a Brownian motion and X^'^^ satisfies the stochastic differential 
equation 

dx\"^ = dW^"^ + v^^\x['''^)dt (5.1) 

where v'^'^\x) = av{ax). This makes it clear that X'^'^^ —>■ as cr 
weakly on C([0, 00), M"^). Recall from Section [U the fact that X^'^'> X 
as cr — > 00, in the same sense, where X is a Brownian motion with 
diffusivity a. 
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Take a = r. Then 

so 

Hence the hmit 

7(wM,i^) := hni l^f 

exists and equals r^~''-j{v,rK). The weak hmits for X^'') as r — or 
r — > oo suggest the following result, which however requires further 
argument because the asymptotic growth rate of the sausage is not a 
continuous function on C([0, oo), M''). Write ca.p{K) for the Newtonian 
capacity of K and capa(if) for the capacity of K with respect to the 
diffusivity matrix d. Thus 

capg(/v) = Vdet d cap{d~^^^K). 

Theorem 5.1 We have 

lim r^~''7(w, rK) = lim 7(1/'') , K) = 7(0, A') = cap(A:) 

and 

lim r'^~''-'-f{v,rK) = lim '-f{v^''\K) =ca.p^{K). 

Proo/ Fix T € (0,oo) and write for the interval ((j - 1)T, jT]. 
Consider for 1 ^ j ^ k the function Fjjt on C([0, 00), R"*) defined by 

Then Fj^k is continuous, so 
lim E(F,- = E(F,- lim E(F,- ..(xW)) = E(i^,, 

r— >0 r— >oo 

Choose Xq so that 7r(Xo) is uniformly distributed. Then by stationarity 
E(Fj-fe(X(''))) = E(Ffc_j(X(''))) where Fj = Fij+i- Fix r and write 
Sf{X^'-'>) = Si. Note that 

5'(0,nT] = 'S'/(i) U • • • U S'/(„). 

So, by inclusion-exclusion, we obtain 

nE(Fo(xM))-^(n-j)E(F,(xM)) ==: E(|5(o.„t] l/T) ^ nE(Fo(xM)). 
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Divide by n and let n — > oo to obtain 

oo 

E(Fo(xW)) - ^E(F_,-(XW)) ^ 7(«W,if) ^ E(Fo(xW)). 

Fix 9 e (1,3/2) and p G (3, oo) with p^^ + q^^ = 1- By Lemmas and 
14.21 there is a constant C(p, u, K) < oo such that, for all r and j, 

E(F,(xM)) = E(|5,(i)n^,(,+i)|/r) 

E ( 1 1^,(1) |/rn'/> (5,(1) n ^ 0)'/' 2C(j - ly"/'^. 

Given e > 0, we can choose J(p, K) < oo so that 

oo oo 

^ E(Fj(xW))s^2C^j-'^/29^^^ (5 2) 

j=j+i j=j 

Wc follow from this point the case r — > oo. The argument for the other 
limit is the same. Let r — > oo to obtain 
,/ 

E(Fo(X)) - VE(F,-(X)) - e ^ liminf 7(i;W, A') 

s$ limsup7(z;(''\ AT) < E(i^o(-^))- (5.3) 

r— >-oo 

It is known that 

Km E(^l,(^)) = lim E(|5f (X)|/T) = cap,(/0. 

T— s-oo T— s-oo 

See Q for the case a ~ I. The general case follows by a scaling trans- 
formation. Note that, for j ^ 1, 

i+i 

\S(0,T]i^) ^ S{JT,U+i)T]i-^)\ + l'5'(^',(i+l)T](-^)l ^ 51 l'^ai-l)T,iT](-^)|- 

i=l 

Take expectation, divide by T and let T — > oo to obtain 

U + l)cap,(/^) +limsupE|5(^^j,,(X) n (X)|/r 

T-i>oo 

< (j + l)cap5(A0 
which says exactly that 

lim E(F,(X)) = 0. 
Hence the desired limit follows on letting T — > oo in (|5.3p . □ 
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6 Homogenization of the advection-diffusion 
equation in a perforated domain 

Our main results are analogues to the homogenization statements (j2.3p . 
(|2.4p for advection-diffusion in a perforated domain. Recall that v is 
a Z'^-periodic, Lipschitz, mean-zero, divergence-free vector field on K'*, 
and K is a compact subset of M.'^. The domain C R'' is constructed 
by removing all the sets K + p, where p runs over the set P of atoms of 
a Poisson random measure /i on M.'^ of intensity p = ct~^. Write 

a = a{v,K), \ ^ X{v,K) ^ j{v,k). 

Theorem 6.1 Let f e L'^(R'^) and a G (0,oo) be given. Denote by u 
the solutioi^ to the Cauchy problem for 

C = iA + w(x).V 

in [0,cxj) X D with initial data f(-/a), and with Dirichlet conditions on 
the boundary of D. Set u^'^\t^x) ~ u{a^t,ax). Then 

E \u'^''\t,x) -u{t,x)\^ dx ^ 

as <T — > oo, for all t ^ 0, where u is the solution to the Cauchy problem 
for i div aV — A in [0, oo) x R'' with initial data f . 

Proof Replace t by cr^t, x by ax and / by f{-/cr) in (|1.2[) and (|1.3p to 
obtain 

E (u^"\t^x)) =E,, (/(x('^))cxp{-p|5f.,(X)|}) 

and 

E(«(-)(t,x)2) =E.. [f{x["'^)f{Y,^^^)eM-p\sX{X)USX{Y)\}) 

where the subscript ax specifies the starting point of X and where Y 
is an independent copy of X. We omit from now on the superscript K. 



^ We extend u to a function on [0, oo) X by setting u[t, x) = for any x ^ D. 
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TheiH 

¥.[\u^^\t,x)-u{t,x)\'^^ 

+ {E,,(/(x(^'))-E,(/(XO)}e-^*)' 

^ E,, E,, n 

+ 2|4"^(t,a;)-uo(i,a:)P (6.1) 

where u^'^ and mq denote the correspondmg solutions to the Cauchy 
problem with initial data / in the whole space, and we used Cauchy- 
Schwarz and (a + bf 2a^ + 25^ and le"" - e^''^ 1^ - a| to obtain 
the inequality. Now 

/ ¥.,Jf{x\''^f) dx^ ( \f{x)\^dx<^ 

because dx is stationary for X'-'^^ and, by (|2.3p . as cr — J' oo 

/ |4''^(t,a;) - ■tio(t,x)pdx -!■ 0. 

So, using Lemma H^l on integrating (|6.ip over R'^ and letting cr — oo, we 
conclude that the right-hand side tends to 0, proving the theorem. □ 

Theorem 6.2 Let a; g E'' and a G (0, oo) be given. Let X be an 
C-diffusion in K'^ starting from ax and set 

T = \ui{t -^Q : Xt(^ K + P}. 

Set X^"^^ = a^^X„2^ and T^'^^ = cr^^T. Write X for a Brownian motion 
in K'' with diffusivity a starting from x, and write T for an exponential 
random variable of parameter X, independent of X . Then, as a oo, 

(X('"),T('")) ^ {X,T), weakly on C([0, oo), R'^) x [0,oo). 

Proof Write 5*4 for the {X, i5r)-sausage. Fix a bounded continuous func- 
tion F on C([0,oo),K^) and fix t > 0. Then 

E (f(X('^))1{^(„,>,j) = E (f(x('^)) cxp{-p|5,2,|}) 

3 This is an instance of the formula E(|X - ap) = var(X) -|- {E(X) - a)^. 
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and 
so 

^ \\F\\ooE^MSaH\ ~ Ai| + |E(F(X('^)) - E(F{X))\e-'^K 

On letting a — > oo, the first term on tlie riglit tends to by Lemma 14.31 
and the second term tends to by (|2.4p , so the left hand side also tends 
to 0, proving the theorem. □ 



7 The case of diffusivity e^I 

In this section and the next we fix e S (0, oo) and consider the more 
general case of the operator 

The following statements follow from the corresponding statements ab- 
ove for the case e = 1 by scaling. Fix x G M'' and let X be an /^-diffusion 
in M.'^ starting from x. Then 

\S^{X)\/t 7(e, V, K), P^-almost surely 

as < — ?> oo, where 

7(£, V, K) = £^7(e"^w, K). 
Moreover, setting v^'^\x) = rv{rx), as above, we have 

^{£,v^'^\K) cap£2/(iir) = £^cap(i4r), as r 

and 

7(£, iC) ^ capa(^_„)(i^), as r ^ oo 

where 

a(£, v) = £^a(£~^t;). 

Fix <T €E (0, oo) and suppose now that X starts at ax. Define as above 
T = inf{t Q : Xt K + P} mid write X^"^ ^ <7^^X^2t and T'^"^ 
a^^T. Then, as cr — ?> oo, 

(X,r), weakly on C([0,oo),M'') x [0,oo) 
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where X is a Brownian motion in of difFusivity a(e, v) starting from 
X, and where T is an exponential random variable independent of X, of 
parameter A(e, w, K) = 7(e, v, K). 



8 Monte Carlo computation of the asymptotic 

growth rate 

Let X be as in the preceding section. Fix T € (0, oo). The foUowing 
algorithm may be used to estimate numerically the volume of the {X, K)- 
sausage St = S!^{X). The algorithm is determined by the choice of three 
parameters N , m, J e N. 

• Step 1: Compute an Euler-Maruyama solution {XnAt : n = 0, 1, . . . , 

N) to the stochastic differential equation 



dXt = edWt + v{Xt) dt, Xq = x 

up to the final time T = NAt (Figure [57Tk ). 
• Step 2: Calculate 



(8.1) 



Rk ~ max \y 

yeK. l^k^d 



IX 



l<n<N. l<k<a 



nAt 



We approximate St by 5^^' = Uo^n^jv (K + XnAt)- Note that s!^^^ 
is contained in the cube with side- length L = 2{Rk + Rx.t) centred 
at X (Figure IHUb) ■ 




(a) Step 1 (b) Step 2 

Figure 8.1 First two steps of the algorithm 



• Step 3: Subdivide the cube of side-length L centred at x into 2*^ sub- 
cubes of side-length L/2 and check which of them have non-empty 
intersection with s!^\ Discard any sub-cubes with empty intersection. 
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Repeat the division and discarding procedure in each of the remaining 
sub-cubes [Figure \8.}3j) iteratively to obtain I sub-cubes of side-length 
centred at yi, . . . , yi say, whose union contains s!^\ 




(c) m = 3, I = 40 (d) m = 4, / = 111 



Figure 8.2 x — y projection of the path and the sub-cubes for different 
values of m 

• Step 4: Generate uniform random variables Ui, Uj in [—1/2, 
1/2]'' and estimate V — \St\ by 

1=1 j = l \ n=0 / 

where Mi,j) = {ly. + 2-™LJ7j + XnAt}- 

The algorithm was tested in the case d = 3. We took e — 0.25 and 
took V to be the Taylor-Green vector field in the first two co-ordinate 
directions; thus 

v{x) = (— sinxi C0SX2, cosxi sinx2, 0)"^. 

We applied the algorithm to X^'^\ which has drift vector field v''^^x) = 
rv{rx), for a range of choices of r e (0,oo). We took K to be the Eu- 
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clidean unit ball B and computed \S^{X^''^)\ for T = 10", using para- 



meter values TV 



10^ 



4 and J = 10^. The numerical method 
used to solve (|8.ip was taken from [ij. The values |S'|(X(''))|/T were 



taken as estimates of the asymptotic volume growth rate 7(0.25, u^'"'', B). 
These are displayed in Figure [ 




Figure 8.3 Growth rate of the sausage for different values of r 

In Section [7] we stated the following theoretical limit, deduced from 
Theorem 15.11 



lim 7(£, v^''\B) = cap(B) = 2TTe'^ = 0.3927. 



(8.2) 



This is consistent with the computed values of 7(0.25, v^^\ B) when r is 
small. 

It is known [13] that a(e, w) has the form 




a(e, v) = 

for some a = Q;(e,w), which can be computed using Monte Carlo sim- 



ulations. In 



16| . this was carried out for e = 0.25 up to a final time 



T = lOf, using a time step At = 10 ^, again using a numerical method 
from PJ to solve The value a{0.25,v) = 0.0942 was obtained as 

the sample average over 10* realizations of |X^p/T. Using this value, 
we simulated X and used the volume algorithm to compute \S^ {X)\/T 
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as an approximation to capa(o.25,ii) (^): obtaining the value 1.4587. We 
showed theoretically that 

lim 7(e,«W,i?)=cap3(,_„)(B). 

The computed value for capg^Q 25 v){-^) consistent with the computed 
values of 7(0.25, B) for large r. 

References 

[1] Aronson, D. G. 1967. Bounds for the fundamental solution of a parabolic 

equation. Bull. Amer. Math. Soc, 73, 890-896. 
[2] Bensoussan, A., Lions, J.-L., and Papanicolaou, G. 1978. Asymptotic 

Analysis for Periodic Structures. Stud. Math. Appl., vol. 5. Amsterdam; 

North-Holland. 

[3] Eisele, T., and Lang, R. 1987. Asymptotics for the Wiener sausage with 

drift. Probab. Theory Related Fields, 74(1), 125-140. 
[4] Fannjiang, A., and Papanicolaou, G. 1994. Convection enhanced diffusion 

for periodic flows. SIAM J. Appl. Math., 54(2), 333-408. 
[5] Hoang, V. H. 2000. Singularly perturbed Dirichlet problems in randomly 

perforated domains. Comm. Partial Differential Equations, 25(1-2), 355- 

375. 

[6] Jikov, V. v., Kozlov, S. M., and Olemik, O. A. 1994. Homogenization of 

Differential Operators and Integral Functionals. Berlin: Springer- Verlag. 

Translated from the Russian by G. A. Yosifian [G. A. losif'yan]. 
[7] Kac, M. 1974. Probabilistic methods in some problems of scattering 

theory. Rocky Mountain J. Math., 4, 511-537. Notes by Macki, J., and 

Hersh, R. Papers arising from a Conference on Stochastic Differential 

Equations (Univ. Alberta, Edmonton, Alta., 1972). 
[8] Kingman, J. F. C. 1968. The ergodic theory of subadditive stochastic 

processes. J. Roy. Statist. Soc. Ser. B, 30, 499-510. 
[9] Le Gall, J.-F. 1986. Sur la saucisse de Wiener et les points multiples du 

mouvement brownien. Ann. Probab., 14(4), 1219-1244. 
[10] Majda, A. J., and Kramer, P. R. 1999. Simplified models for turbulent 

diffusion: theory, numerical modelling and physical phenomena. Physics 

Reports, 314, 237-574. 
[11] Norris, J. R. 1997. Long-time behaviour of heat flow: global estimates 

and exact asymptotics. Arch. Rational Mech. Anal., 140(2), 161-195. 
[12] Papanicolaou, G. C, and Varadhan, S. R. S. 1980. Diffusion in regions 

with many small holes. Pages 190-206 of: Stochastic Differential Systems 

(Proc. IFIP-WG 7/1 Working Conf, Vilnius, 1978). Lecture Notes in 

Control and Inform. Sci., vol. 25. Berlin: Springer- Verlag. 
[13] Papanicolaou, G. C, and Varadhan, S. R. S. 1981. Boundary value 

problems with rapidly oscillating random coefficients. Pages 835-873 of: 



Homogenization for advection- diffusion in a perforated domain 19 

Random Fields, Vol. I, II (Esztergom, 1979). CoUoq. Math. Soc. Janos 

Bolyai, vol. 27. Amsterdam: North-Holland. 
[14] Pavliotis, G. A., Stuart, A. M., and Zygalakis, K. C. 2009. Calculating 

Effective Diffusivities in the Limit of Vanishing Molecular Diffusion. J. 

Comp. Phys, 4(228), 1030-1055. 
[15] Ranch, J., and Taylor, M. 1975. Potential and scattering theory on wildly 

perturbed domains. J. Fund. Anal, 18, 27-59. 
[16] Zygalakis, K. C. 2008. Effective Diffusive Behavior For Passive Tracers 

and Inertial Particles: Homogenization and Numerical Algorithms. Ph.D. 

thesis. University of Warwick. 



